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ABSTRACT 



Aims. The two dimensional structure of hot gas in galaxy clusters contains information about the hydrodynamical state of the cluster, 
which can be used to understand the origin of scatter in the thermodynamical properties of the gas, and to improve the use of clusters 
to probe cosmology. 

Methods. Using a set of hydrodynamical simulations, we provide a comparison between various maps currently employed in the X- 
ray analysis of merging clusters and those cluster maps anticipated from forthcoming observations of the thermal Sunyaev-Zel'dovich 
effect. 

Results. We show the following: 1) an X-ray pseudo-pressure, defined as square root of the soft band X-ray image times the tempera- 
ture map is a good proxy for the SZ map; 2) we find that dumpiness is the main reason for deviation between X-ray pseudo-pressure 
and SZ maps; 3) the level of dumpiness can be well characterized by X-ray pseudo-entropy maps. 4) We describe the frequency of 
deviation in various maps of clusters as a function of the amplitude of the deviation. This enables both a comparison to observations 
and a comparison to effects of introduction of complex physical processes into simulation. 

Key words. Galaxies: clusters: general - Cosmology: theory - Cosmology: observations - X-rays: galaxies: clusters 



1. Introduction 

The X-ray emission from clusters of galaxies has been conven- 
tionally characterized in terms of an X-ray surface brightness 
image and a temperature map. Historically, due to the success 
of the Einstein (Gursky et al. 1972) and ROSAT (Boehringer et 
al. 2000, 2001) missions, X-ray imaging data are available for a 
wide range of clusters, both in mass and redshift, while cluster 
searches based on the level of X-ray emission have become one 
of the primary cosmological tools (for a review see Borgani & 
Guzzo2001). 

The availability of cluster temperatures, on the other hand, 
has long been limited to single emission-weighted values, as 
provided by Ginga, Einstein and EXOSAT (David et al. 1993; 
Edge et al. 1992), and although the first temperature maps 
were made with ROSAT (Briel & Henry 1994) and ASCA 
(Markevitch 1996, 1999) data, they were limited in either spec- 
tral range (ROSAT) or spatial resolution (ASCA). With the ad- 
vent of Chandra and XMM-Newton, temperature maps have be- 
come widely available, revealing an unprecedented amount of 
detail (Vikhlinin et al. 2001; Mazzotta et al 2002; Briel et al. 
2004; Finoguenov et al. 2004; Henry et al. 2004). However, nei- 
ther temperature maps nor X-ray images are the primary char- 
acteristics of the gas, and are instead results of underlying pres- 
sure and entropy distribution, that trace the dark matter poten- 
tial, reflects the thermodynamical history of the gas and the on- 
going magneto-hydrodynamical processes. However, projection 
effects complicate this study and in practice X-ray pseudo pres- 
sure and entropy maps have been introduced (Markevitch et al. 



1996; Finoguenov et al. 2004, 2005, 2006a,b, 2007; Briel et al. 
2004; Henry et al. 2004), addressing an immediate goal of dis- 
criminating between the various origins of fluctuations, e.g. to 
differentiate between shocks and gas displacement, also known 
as cold fronts (Vikhlinin et al. 2001). The way these maps work 
is that at the position of the cold front the contrast in the entropy 
map is increased, while the contrast in the pressure map is re- 
moved, compared to either images or temperature maps, while 
on the positions of the weak shocks, the contrast in the entropy 
maps is suppressed, while the contrast on the pressure map is 
enhanced. Some astonishing discoveries using these maps were 
made using very deep XMM observations of M87, recovering 
the details of AGN feedback (Simionescu et al. 2007) and the 
underlying triaxiality of the dark matter distribution. 

These pressure and entropy distributions are called pseudo 
maps to reflect the fact that both pressure and entropy are lo- 
cal quantities, while the projected emission weighted tempera- 
ture and X-ray surface brightness are averaged along the line 
of sight. The effect of projection in reducing the fluctuation has 
been found to be a function of the spatial scale of fluctuation 
(Schueckeret al. 2004), with small-scale fluctuations being more 
strongly suppressed. Since all cluster maps suffer from projec- 
tion effects, it is important to understand the link between vari- 
ous methods of studying clusters, such as thermal SZ, weak lens- 
ing and X-ray. 

In order to improve our understanding of cluster physics, we 
employ a set of hydrodynamical simulations to study how well 
various cluster maps compare. We paid a particular attention on 
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X-ray image 




Fig. 1. Example of the cluster maps in one of the simulation runs. From the top left to the bottom right we show the dark matter 
column density, gas column density, X-ray images in the 0.5-2 keV band, X-ray temperature map, pseudo entropy map, pseudo 
pressure map, Compton Y map, dumpiness map, ratio of Compton Y to X-ray image. The grey scale is set to depict the central 
structure in black in all the maps. 



possible use of cluster maps for diagnostics not readily avail- 
able to the observer. In §2, we describe the simulation sample. In 
§3 we look for correlation among the maps, discuss the general 
trends and look for correlations of deviations from the general 
trends. We also present the statistical expectation for deviations 
to occur as a function of their strength. We summarize our find- 
ings in §4. 



2. Mock Observations 



We use a catalog of 68 galaxy clusters evolved with the P3MSPH 
code (Evrard 1988, 1990), incorporating preheating of the ICM 
to mimic the effects of galaxy feedback. The level of entropy 
introduced into the initial conditions is tuned to 106 keV cm 2 
in order to match the observed scalings of luminosity and ICM 
mass with temperature (Bialek et al. 2001). No additional AGN 
feedback and radiative cooling has been included into the simu- 
lations. 
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Fig. 2. The distribution of correlation coefficients between various simulated maps for each of the 68 simulated clusters, calculated 
using EquationQ] Left panel represents correlations between dark matter column density and each of the pressure (solid), Compton 
Y (dashed) and square root of the X-ray surface brightness 7^ 5 (dotted). Right panel represents correlations between Compton Y 
and each of the pressure (solid) and I^ 5 (dotted). The dashed line shows a correlation between the ratio of Compton Y and X-ray 
surface brightness with pseudo entropy. All comparisons are carried out within 1/4 of r2oo- 



The simulations were produced using a multi-step procedure 
outlined in (Bialek et al. 2001). A single cluster from this dataset 
was found to exhibit a post-merger cold front/features driven by 
sub-halo sloshing, as presented in Bialek, Evrard & Mohr (2002) 
and further detailed in Pawl et al. (2005). The simulations evolve 
isolated cut-outs of a larger N-body realization, with regions 
chosen such that final halo masses are sampled with roughly 
equal likelihood across the mass range M200 = 0.2 - 2 x 10 15 M Q . 
The runs are of moderate resolution, with gravitational soften- 
ings between 50 and 150 kpc, and dark matter particle masses 
of 0.4 - 4 x 10 10 M o . While Lewis etal (2000) caution that large 
softening can affect the overall cluster potential, the comparison 
study of Frenk etal ( 1 999) demonstrates that moderate resolution 
runs offer good convergence in dark matter and gas properties on 
scales larger than roughly twice the hydrodynamic smoothing 
scale. 

A typical halo is resolved by ~ 50, 000 particles. With pre- 
heating of 106 keV cirT 2 , this is sufficient to resolve low-order 
measures of ICM properties, such as temperature and luminos- 
ity, and morphological features driven by mergers with mass 
ratios of order 0.1 and larger. The underlying cosmology is a 
flat, concordance model with Q m = 0.3, Q.a = 0.7, Q./, = 0.03, 
<t 8 = 1.0, and h = 0.7, where the Hubble constant is defined as 
lOO/z km s Mpc~' and erg is the power spectrum normalization 
on 8/T 1 Mpc scales. Details of the full ensemble are presented 
on the project web-page 1 . 

In this paper we only used the final (z=0) configuration 
of 68 model clusters and employ 3 orthogonal projections of 
each model as independent data, giving a total of 204 separate 
datasets. Every map has 256 pixels on a side, spanning a re- 
gion twice the virial radius of the system. By construction the 
maps are designed to study the appearance of individual clus- 
ters and not the projected effects in cluster surveys (e.g. White 
2003). From an initial set of X-ray surface brightness (I) and 
temperature maps (T), we have constructed the X-ray pseudo 
pressure and entropy maps (hereafter P and S, respectively) as 



p = I 1 1 2 * T and S = T/I 1 ^ 3 and compared those to the column 
density of the both ionized gas (G) and dark matter (DM), as well 
as the predicted Compton Y parameter of the thermal Sunyaev- 
Zel'dovich (SZ; Sunyaev & Zel'dovich 1980) effect (Y). We de- 
fine the dumpiness as C — G/ yl, which differs from its exact 
definition by the square root of the projection length. In Fig. Q] 
we show the simulated cluster maps obtained for the gl run, 1st 
projection. 

The dark matter maps in simulations appear much clumpier 
compared to the gas. This is likely due to erosion of small-scale 
structure in the gas by setting the entropy floor. A careful com- 
parison between dark matter subhalos and the hot gas has been 
made in Powell et al. (2009). Using the case of one halo sim- 
ulated with and without cooling, they show that occupation of 
small subhalos by hot gas is efficiently suppressed under both 
treatments. Dark matter maps of cluster are thus likely to ap- 
pear clumpier than the X-ray gas irrespective of the detailed gas 
physics assumptions. Near equilibrium, the gas will trace the 
equipotential surfaces, which tend to be smoother and rounder 
than the density isosurfaces. In order to remove the small-scale 
power in the comparison, we have convolved the dark matter 
maps with the Gaussian. The members of the simulated sample 
range in spectral temperature from 1.5 keV to about 8 keV, with 
cluster masses M200 ranging from (0.015-2.4) xlO 15 M . The de- 
rived images, along with associated parameters describing the 
simulations are publicly available as part of VCE, the Virtual 
Cluster Exploratory . 



3. Results 

For each simulated cluster we can quantify the correlation be- 
tween the different maps within a fixed fraction of the virial ra- 



1 http://vce.physics.lsa.umich.edu 
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dius. We use the following expression for the correlation coeffi- 
cient between images X and Y for a given cluster 

r= _ = =, (1) 

where i sums over all the pixels in the image. We search 
for correlations between all primary maps. The distribution of 
correlation coefficients across the ensemble of simulated clusters 
is shown in Fig. [2] The strongest correlation occurs between P 
and Y. The dark matter is equally well correlated with either Y 
or P, and somewhat less well with I. 

We have found strong correlations between X-ray pseudo- 
entropy and the ratio of Y to I; which we attribute to the fact that 
the entropy can also be defined as Y/I 5 ^ 6 . This can useful in case 
where only X-ray imaging is available as well as good SZ-data, 
but no temperature map, as should be the case for high-z clusters. 
Historically, the ratio of Y-to-I was proposed for measurement of 
the Hubble constant (see e.g. Bonamente et al. 2007). Our results 
provide an insight on the origin of scatter in this measurement, 
and our proposal to reduce the scatter in this comparison consists 
of the use of pseudo pressure maps instead of X-ray images and 
also estimating the expected deviation based on the amplitude of 
the scatter seen in the pseudo entropy maps. 

3.1. Estimated influence of the general trend 

In this section we estimate analytically the effect of projection 
on the derived correlation coefficient, accounting therefore for 
differences in the radial trends between the maps and compare 
these predictions to the observed correlation coefficients. 
Consider the scaling with radius of different properties. 




A/3 in (l+( r / r js)-(o.e3 + m 



Fig. 3. Variation in the correlation coefficient as a function of 
change in the radial slope. The changes are calculated in respect 
to the Compton Y two-dimensional distribution with the profile 
of (1 + (r/r c ) 2 )~ <0 - 62+A/3) . The points show the actually measured 
correlation strength and the range enclosing 68% of the data. 
The dotted, dashed and solid lines represent a correlation be- 
tween Compton Y and each of the pseudo pressure, the dark mat- 
ter column density and the square root of X-ray surface bright- 
ness (/^ 5 ), respectively. While the pseudo pressure and Compton 
Y parameter exhibit a correlation consistent with the predicted 
trend, the correlation with both dark matter and the X-ray image 
is lower than the expectation based merely on differences in the 
profiles. 
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where n is the gas density, is the radial slope for which 
we assumed a typical cluster value. The entropy profile is as- 
sumed, based on the results of observations (e.g. Finoguenov et 
al. 2005). The DM is the projected column density of dark mat- 
ter, for which the distribution was taken to resemble closely the 
NFW profile (Navarro, Frenk & White 1997). Other parameters 
can be derived based on these assumptions. Addition of the con- 
stant 0.5 is a result of the Abell integral. 

In order to illustrate the effect of differences in the ra- 
dial slope on the observed correlation, in Fig.. [3] we calculated 
an expected reduction in the correlation coefficient, assuming 
(1 + (r/r c ) 2 )~ (a62+A/3) , varying A/3 to cover the range of all the ob- 
served clusters. The main conclusion of this test is that while the 
typical correlation coefficients between Compton Y vs. pseudo 
pressure are reproduced, the observed correlations between sim- 
ulated X-ray image and Compton Y as well as all the correla- 
tion of Compton Y with dark matter column density are much 



weaker. This implies both that the gas is only an approximate 
tracer of the dark matter and that hydrodynamical processes in- 
volving gas displacements (cold fronts) are playing a dominant 
role in defining the soft X-ray band images of the cluster, in ac- 
cordance with recent Chandra and XMM results (Vikhlinin et al. 
2001, Mazzotta et al. 2002; Briel et al. 2004) and other numeri- 
cal simulations (Powell et al. 2009). 

3.2. Definition of substructure as deviation from the mean 
trend 

In order to see the correlation in the substructure on the im- 
ages, we have corrected for the radial variation in all maps. We 
perform this correction by calculating an average map for each 
quantity, across the ensemble of simulated clusters, which cap- 
tures the broad systematic trend. This average map is then used 
to normalize each individual cluster map, to reveal the devia- 
tions about the global trend caused by interesting cluster physics. 
The correlation between these corrected maps then more cleanly 
probes the relation between these deviations. 

To eliminate bias from anomalous clusters, we identify and 
exclude outliers from the calculation of the average map using 
the following iterative scheme. Firstly, we take the logarithm of 
all the pixel values in the map for every cluster and calculate 
the mean of these logged images. We then compute the correla- 
tion coefficient between this mean logged image and the individ- 
ual logged map for each cluster. Any clusters with a correlation 
below 0.85 are identified as outliers and the mean logged im- 



A. Finoguenov et al.: Hydrodynamic simulations of correlation and scatter in galaxy cluster maps. 



5 



age is recalculated with them excluded; the process is repeated 
until no further clusters are excluded. This process is found to 
change the resultant map by only ~5%. The distribution of cor- 
relation coefficients for the corrected maps is plotted in Fig. [4] 
and a comparison with Fig. [2] shows the much broader spread 
of values obtained when the mean trend is removed. Intuitively, 
the level of correlation between the gas-related maps and dark 
matter maps should be quite sensitive to the physics introduced 
in simulations, such as cooling, conduction and viscosity. Thus 
a comparison of predictions of various simulations to the ob- 
servational data shall be rewarding. The correlation between the 
gas-tracing maps has also decreased. For example a correlation 
between the pseudo-pressure to Compton Y maps dropped from 
0.98 to 0.95. In order to understand the origin of this decline 
we move to our next test - a direct comparison of the levels of 
fluctuations of each cluster. 

3.3. Scatter 

In addition to the similarity of the maps as traced by the corre- 
lation estimators, we can quantify the link between the maps by 
examining the degree of deviations in the map ratios. It is also 
advantageous to compare the deviations seen in the maps to devi- 
ations from the average trend, as well as their typical amplitudes. 
Our primary results are outlined in Fig. [5] 

When comparing the rms level for pseudo pressure and 
Compton Y, one sees that although the two scale very well, the 
pseudo pressure depicts the large deviations more strongly. This 
explains why there is a reduction in the correlation strength be- 
tween taking the full map or only comparing the deviations from 
the mean trend. A comparison between the pressure and entropy 
deviations shows that although the pressure deviations are typi- 
cally larger, if all the deviations would be due to the shocks, the 
expected level of entropy deviations would be limited to 0.05. 
The observed entropy fluctuations are typically 10 times larger 
and therefore require different explanation, such as survival of 
low-entropy gas during the cluster assembly. A similar range of 
the fluctuations is reported for the REFLEX-DXL cluster sample 
(Finoguenov et al. 2005); in Poole et al. (2007, 2008) the ampli- 
tude of fluctuations in the maps is further linked to the merger 
state of the cluster, which has also been suggested observation- 
ally by Zhang et al. (2009). 

Fig- reveals a good correlation between the dispersion in 
the entropy and dumpiness. The very presence of such a corre- 
lation implies that the dumpiness of the gas is dominated by the 
large-scale gas displacements, which are traced by the entropy 
map. Previously, we have shown that the entropy deviations are 
nearly cospatial with deviations in dumpiness, as shown by their 
reasonable correlation shown in Fig. [4] It implies that to reduce 
the effect of dumpiness (which for example introduce a scatter 
in the measurement of Hubble constant based on the X-ray and 
SZE comparison), one should excise the zones most deviating in 
the entropy maps. 

This suggestion can be further strengthened when consider- 
ing the rms of the P-to-Y ratio. First, we note that the rms in the 
individual pressure maps are much higher than the rms of the 
P-to-Y ratio. The latter has a scatter lower by a factor of 5, typi- 
cally at the 10% level. To prove that dumpiness is the reason for 
these residual fluctuations, we first note that a ratio between P 
and Y could be written as rx y ^ = where C is dumpiness. 
So we construct the additional map, C p = -^y^ to examine the 
fluctuations associated with averaging in Compton Y maps. The 
scatter in P/Y is a factor of two larger than in C p , which proves 



that the dumpiness is the dominant effect in producing devia- 
tions between pressure inferred from X-ray observations and SZ 
signal. 

To estimate the importance of certain level of deviations, we 
next perform the analysis of the frequency of their occurrence. 
Fig.|6]displays the histogram of the rms distribution for the ma- 
jor parameters of this study. Already in Finoguenov et al. (2005, 
2007) we have compared the observations to these predictions, 
revealing deviations in pressure fluctuations observed on the 
mass scale of galaxy groups. The first distribution of dark matter 
substructure has been presented in Smith et al. (2005), mani- 
festing that comparison with various datasets will not be limited 
to X-ray observations. As evident from Fig. |6j the Compton Y 
has the smallest deviations compared to either dark matter, X- 
ray pressure or the square root of the X-ray image. The latter 
three exhibit similar deviations, with the exception that for the 
dark matter, deviations at the 25% level are seen much more fre- 
quently. As already mentioned above, the deviations in P/Y are 
intermediate between the temperature and C p . Fluctuations in 
the dumpiness itself are much larger and the frequency of their 
appearance is similar to that of entropy. 

4. Conclusions 

We have used a set of hydrodynamic al galaxy cluster simulations 
to investigate the two-dimensional comparison between maps of 
different quantities which are increasingly becoming available 
to the observer. We test the idea of using these maps for hydro- 
dynamical analysis and for achieving better accuracy in using 
clusters for cosmological studies. We conclude that the X-ray 
pseudo pressure map is a very good proxy for the Compton Y 
image obtained from SZ observations. It could be used to simu- 
late the observational appearance of the clusters on the SZ sky, 
their scatter from the averaged trend and, most importantly, pro- 
vide the lowest scatter comparison between spatially resolved 
SZ and X-rays measurements. 

We have also shown that the pseudo entropy maps can be 
used as a proxy of dumpiness of the gas, which cannot otherwise 
be observed directly. We demonstrate that the scatter in the ratio 
of the X-ray pressure estimate to the Compton Y is driven by 
dumpiness and not by shocks, so elimination of clumpy zones, 
e.g. marked as deviant on the entropy maps, is the most criti- 
cal step towards achieving a consistent picture of clusters with 
various methods. With the prescriptions discussed in this work, 
we anticipate a reduction of scatter in the SZE-to-X-ray com- 
parisons from 30% to 5%, through the use of X-ray pressure 
restoration and excising strongly deviant zones in the entropy 
maps. 

We have shown that the structure in dark matter maps are 
poorly correlated with the structure in any gas-based maps. 
Therefore the power of using the gas-maps to trace the dark mat- 
ter substructure is quite limited. A degree to which this correla- 
tion will be observed in the real data is sensitive to the detailed 
physics of gas disruption and as such should allow to constrain 
the role of cooling, conduction and viscosity through matching 
the observations to sets of hydrodynamical simulations, as ad- 
vocated in Dolag et al. (2004). We have also introduced a dif- 
ferent test, which describes the frequency of the observed rms in 
the maps. The has already enabled a comparison to observations 
(Finoguenov et al. 2005; 2007), revealing that while the scatter 
in both entropy and pressure maps in clusters agree well with 
the prediction of these simulations, the amount of scatter seen in 
the pressure maps of galaxy groups is much larger, indicating a 
larger role of internal feedback processes there. 
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Fig. 4. The distribution of the pixel -by-pixel correlation coefficient (Equation Q]) between the maps, after correcting for the mean 
trend in each map. Left panel represents correlations between dark matter column density and each of the pressure (solid), Compton 
Y (dashed) and square root of the X-ray surface brightness /J 5 (dotted). Middle panel represents correlations between Compton Y 
and pressure (solid) and /£- 5 (dotted). Right panel represents correlations between the entropy and the dumpiness (solid line) and 
correlations between the ratio of Compton Y and X-ray surface brightness with entropy (dashed). All comparisons are carried out 
within 1/2 of rjoo- Changing the radii of the comparison to 1/4 r2oo does not change these results. 
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Fig. 5. Correlations between the rms values of various maps with respect to the mean trend. Each point corresponds to a single 
dataset from our ensemble of 68x3 simulated clusters (see Section[2]i. All calculations have been carried out within 0.5 r2oo- The 
dotted lines on all the plots show 1 : 1 correspondence in rms level. 




Fig. 6. Fraction of the sample with rms of studied parameter greater than the jf-axis value. Left panel: shows pressure (grey line), 
square root of X-ray image (solid line), Compton Y (dotted line), dark matter (dashed line). Right panel: shows P2Y (grey), tem- 
perature (solid), dumpiness (dotted), C p (short-dashed), Y2I (long-dashed line), entropy (dot-dashed line). 
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